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Abstract: In this paper we develop a dynamical model of a rotor and the active magnetic bearings 
used to support the rotor. We use this model to develop a stable state feedback control of the magnetic 
bearing system. 

We present the development of a rigid body model of the rotor, utilizing both Rotation Matrices (Euler 
Angles) and Euler Parameters (Quaternions)* In the latter half of the paper we develop a stable state 
feedback control of the actively controlled magnetic bearing to control the rotor position under imbalances. 

The control law developed takes into account the variation of the model with rotational speed. We show 
stability over the whole operating range of speeds for the magnetic bearing system . Simulation results 
are presented to demonstrate the closed loop system performance. We develop the model of the magnetic 
bearing , and present two schemes for the excitation of the poles of the actively controlled magnetic bearing. 

We also present a scheme for averaging multiple sensor measurements and splitting the actuation forces 
amongst redundant actuators. 

1 Introduction 

Several representations of rigid body rotations, including Rotation Matrices , Cayley -Kline parameters , Euler Pa- 
rameters Spinors ( [K.W86] [PP80] [OB79] ) have been developed. Conventional methods of deriving rigid body 
dynamics utilize the Euler angle parametrization of the space of orientations of the rigid body. Such a parametrization 
of 50(3) suffers from coordinate singularities. The singularities are entirely a result of the choice of parametrization. 
A parametrization that is globally nonsingular is the parametrization utilizing Euler parameters (unit quaternions). 
In this paper we develop the dynamical equations describing the rigid body model of a rotor supported by actively 
controlled magnetic bearings, using rotation matrices (parametrized by euler angles) and using euler parameters. 

We begin this paper by giving a brief description of the various mathematical terms and ideas that will be used in 
defining rotation matrices and euler parameters [YCBDB82]. We present some of the properties of rotation matrices 
and quaternions. We derive the dynamical equations of the rotor supported by active magnetic bearings using both 
rotation matrices and quaternions. We present the development of a stable state feedback control law and simulation 
results of the system when controlled by this state feedback control law. Most magnetic bearing systems are comprised 
of redundant sensors and actuators. We present a linear algebraic technique utilizing the method of least squares to 
average multiple measurements and split the actuation forces amongst redundant actuators. 

2 Preliminaries 

The magnetic bearing system utilizes many frames of reference, in which various quantities such as positions, velocities 
and angles are described. In this section we will set out the notation by which we will refer to the various quantities. 
Also we will present some of the definitions and facts necessary for the derivation of the dynamical equations. 
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2.1 Notation 

In this section we will set out the notation by which we will refer to the various quantities. 

• Vectors will be referred to by lower case letters, with an arrow on top. Position vectors will be represented in 
either Cartesian or spherical coordinate systems. Representation of vectors in Cartesian coordinate systems 
will be the components of the vector. Vectors would also be represented as column matrices. The 

components of the column matrix would not contain an arrow. 

Example 2.1 The vector a would be represented in either of the two ways . 
a = CLx $ T a yY + a z 2 


The components of a vector would be written without the arrow on top. 

Matrices and tensors will be referred to by upper case letters. Dimensions and components of the matrix will 
always accompany the notation. Example: A 6 9£ 3x3 

Variables and constants will be denoted by lower case letters, and will be accompanied by a statement concerning 
their dimension. Example: 8 6 S. 

When referring to variables, which are described in a frame of reference, the subscript of the variable will refer 
to the frame of reference. Example: 

We will use three frames of reference primarily. 

1 The first frame of reference is the fixed inertial frame of reference. This frame will be referred to as the 
inertial frame, and the subscript inertial will accompany variables described in this frame of reference. 
The orthonormal coordinate vectors of this reference frame will be denoted as Xi nerOa(> Vir»ertia/» 2; inertial • 

2. The second frame of reference is rigidly attached to the center of mass of the rotor, and moves with the 
rotor. This frame will be referred to as the rotor based reference frame, and the subscript rotor 
will accompany variables described in this frame of reference. The orthonormal coordinate vectors of this 
reference frame will be denoted as Jtrotor , Vrotor, Vrotor- The orientation of the rotor frame is along the 
principal axes of inertia of the rotor. The origin of the rotor frame is at the center of mass of the rotor. 

3. The third frame of reference is attached to the center of the magnetic bearing. We will say more about 
this frame later. The origin of this frame is coincident with the center of the bearing. This frame will be 
referred to as the bearing based reference frame, and the subscript bearing will accompany variables 
described in this reference frame. The orthonormal coordinate vectors of this reference frame will be 
denoted as Afcearinj) Tearing) ^enring- 

Rotation matrices relate the orientation of vectors in one frame relative to another. The convention we employ 
through this report would be as follows. 


L J frame - 1 L J frame-2 

where /jam'll is a rotation matrix that expresses the basis vectors of frame - 2 in terms of the basis vectors 
of /rame-T The rotation matrix can be expressed as a combination of basis vectors of both the frames in 
the following manner. Given that the basis vectors of frames 1 and 2 are represented with the appropriate 
subscript, 


2 ’ A /rame — 1 

y f rame—2 ' 

A f rame-1 

Z frame — 2 * 

^/rame-1 

-2 * ^ frame — 1 

^ frame — 2 

‘ Vyrame- 1 

Z f rame— 2 

' y f rame — 1 

-2 * Z J rame-1 

^ f rame — 2 

• Z f rame— 1 

Z f rame — 2 

‘ ^/rame-1 


• Position vectors of objects will be referred to in the following 

r object 1 


manner. 


' x— frame 

— object object 

^ frame T y — frame 

object 

_ r z — frame 
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The cross product form of a vector is referred to as S(-). That is 
a X b = S(a)b 

0 — a z a v 


S(a) = 


a = 


b = 


a z 

dy 

a x 

a y 

L dz 

b x 


0 


-a 2 

0 


A unit vector in the direction of an object will be represented as 


■y object 


U 


frame 


object 
a x- frame 
object 


U 


y-/ram< 

u ob} ' ct , 

z — f rame J 
-* object 
' frame 


1^‘IU 




/ r obj ect 
= \/[ r *-/ran 


( 4 ) 

( 5 ) 

( 6 ) 


i f object i 2 . [ r °>>jec‘ -12 

' t'y — framei 'l 1 z — frame! 

• By norm, we refer to the Euclidean norm of the vector throughout this report. 

• If the object being referred to is the origin of a coordinate frame, it will be referred to as f/ rame - 2 -origtn 

7 /rame — 1 

To set the ideas clear, consider the following examples. 

Example 2.2 The origin of the rotor coordinate frame, as observed in an inertial reference frame would be 
represented as 


-rotor — origin 

7 inertial 


rotor — origin 
x — inertial 
rotor — origin 
y — inertial 
rotor —on gm 
z — inertial 


(7) 


We tag the word origin when we explicitly refer to the origin of a coordinate frame. 

Homogenous Transforms 

Mappings between points in the Euclidean group 5£’(3) to points in SE( 3) are represented as 4 x 4 matrix 
transformations that map a position and orientation of a frame to another position and orientation. Such 
mappings are termed homogenous transforms , and in coordinates are specified as follows. 

TfcZ-i ■ sm - SE( 3) 


A frame — 1 


A frame— 2 

I/rame-1 

r-pf rame — 2 

^ frame— 2 

& frame — 1 

* frame — 1 

^ f rame— 2 

1 


1 


rpframe—2 
frame — 1 



frame— 2 — origin “ 

* 

nframe— 2 

7 x — frame — 1 
f rame— 2— origin 


**"f rame — 1 

y — frame — l 


[ 0 0 o ] 

f rame— 2 — origin 
_ 1 z — frame — 1 

1 



The inverse of a homogenous transform T J ™™*- 1 represented and given as follows. 


zrpframe—2'\ — \ 
frame — 1 J 


rpframe—1 
* frame — 2 


r pframe — 2iT 
1 frame — 


0 0 0 ’ 


— r fi»/ r «me-2 1 T 
rame — 1 J 


f rame — 2 — origin 
' x — frame — 1 
frame— 2 — origin 
1 y — frame — 1 
frame — 2 — origin 
' z — frame — 1 

l 
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Example 2.3 The homogenous tranformation tranforming coordinates of the origin of the bearing coordinate frame 
as observed in the inertial coordinate frame to the same coordinates as observed from the rotor reference frame would 
be represented as follows. 



( 8 ) 

(9) 


2.2 Definitions 

We will use the properties of linear vector spaces, quaternions and quaternion algebra. In this section we begin by 
defining vector spaces and algebras. We then proceed to state/derive the properties of rotation matrices and euler 
parameters (quaternions and quaternion algebra) [L.A79] [AR88] [AR89] [RA87] [K.W86]. 

Groups: A group is a set X with an internal operation X x X —* ► X, such that 

• the operation is associative 

(xy)z = x(yz) Vz,y,z £ X , 

• there is an element e £ A' called the identity such that 

xe = ex = x V x € X> 

• for each x £ X there is an element of X called the inverse of x (written x" 1 ), such that 

-l -l 

x x = xx — e. 

Usually this group operation is referred to as multiplication. If the operation is commutative then it is referred 
to as addition and the group is called an Abelian Group. 

Ring: A ring is a set X with two internal operations called multiplication and addition, such that 

• X is an abelian group under addition, 

• multiplication is associative and distributive with respect to addition. 

If the group has an element e £ X such that ex = xe = x V x £ X it is called a ring with identity. Also, if 
x £ A' has an inverse, then it is said to be regular. 

Field: A field is a ring with identity, all the elements of which (except zero, the additive identity) are regular. 
Module: A module X over a ring R is an abelian group A" with an external operation, called scalar multiplication, 
such that 

a(x + j/) = otx + ay 

(a 4- fi)x = ctx 4* fix 
( otp)x = a(fix) 
for all a, ft £ 9f? and x, y £ X . 

Algebra: An algebra A is a module over a ring R with identity with an internal operation called multiplication such 
that 

• A is a ring, 

• the external operation (or,z) — ► ax is such that 

a(xy) = (ax)y = x (ay). 

Vector Spaces; A vector or linear space is a module for which the ring of operators is a field. Its elements are called 
vectors. 
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Quaternion: Quaternions can be viewed in many ways. A quaternion is defined as an operator with a scalar q 0 and 
a vector part q, expressed either as a sum of its parts, 

q = l</o *f 

or as a four dimensional vector, 



If go = 0 then the quaternion is called a vector quaternion , and if q = 0 then it is said to be a scalar quaternion. 
In this paper we notate the quaternions by boldface letters. 


2.3 Properties of Rotation Matrix 

The rotation matrix that relates the orientation of one frame relative to another requires the specification of three 
angles, and can be parametrized in a number of ways. We now indicate two commonly utilized parametrizations. 

• Fixed Axis Rotations: Let frame — 1 and frame — 2 be coincident to begin with. 

— Rotate frame — 2 through an angle 9 X about the vector J?/ r ame- 1 . 

- Rotate frame — 2 through an angle 9 y about the vector f/ r am«-i- 

— Rotate frame — 2 through an angle 9 Z about the vector f/ rame _i, 

The resulting rotation matrix, relating the coordinate vectors of frame - 2 to the coordinate vectors of frame — 1 
can be given as 


cos 9 z 

— sin 9 Z 0 


COS 9y 

0 

sin 9 y 


■ 1 

0 

0 

sin 9 Z 

cos 9 X 0 


0 

1 

0 


0 

cos 9 X 

— sin 9 X 

0 

1 


— sin 9 y 

0 

COS 9 y 


. 0 

sin 9 X 

cos 9 X 


Comment 2.1 Note that as the successive rotations are performed about the fixed axes, the rotation matrices are 
premultiplied in the order in which the rotations are performed. 


R, f 


rame—2 
rame — 1 


COS S z COS 9 y 

sin 9 Z cos 9 y 
— sin 9 y 


cos 9 Z sin 9 y sin 9 X — sin 9 Z cos 9 X 
sin 9 z sin 9 y sin 9 X -f cos 9 Z cos 9 X 
cos 9 y sin 9 X 


cos 9 Z sin 0 y cos 9 X + sin 9 Z sin 9 X 
sin 9 z sin 9 y cos 9 X — cos 9 Z sin 9 X 
cos 9 y cos 9 X 


i 

1 

= R l . ram '- 2 l 

" X " 
Y 

L z . 

/rame — 1 

frame — 1 

Z _ 


( 10 ) 


( 11 ) 


• Moving Axis Rotations: Let frame — 1 and frame — 2 be coincident to begin with. 

- Rotate frame — 2 through an angle 9 Z about the vector 2 frame- 2 * 

- Rotate frame — 2 through an angle 9 y about the vector Y/rame- 2- 

- Rotate frame — 2 through an angle 9 X about the vector X f rar Ac- 2 - 


The resulting rotation matrix, relating the coordinate vectors of frame — 2 to the coordinate vectors of frame — 1 
can be given as 


Rf f 


rame. —2 
rame — 1 


cos^ x 

— sin 9 z 

0 ' 


COS 9y 

0 

sin 9 y 


' 1 

0 

0 

sin 9 x 

cos 9 X 

0 


0 

1 

0 


0 

cos 9 X 

— sin 9 X 

0 


1 


— sin 9 y 

0 

COS 9 y 


. 0 

sin 9 X 

cos 9 X 


Comment 2.2 Note that as the rotations are performed about the moving axes, the rotation matrices are post- 
multiplied in the order in which the rotations are performed. 


R 


} 

f 


rame 

rame 


2 

1 


cos 9 Z cos 9 y 
sin 9 z cos 9 y 
— sin 9 y 


cos 9 Z sin 9 y sin 9 X — sin 9 Z cos 9 X 
sin 9 2 sin 9 y sin 9 X + cos 9 Z cos 9 X 
cos 9 y sin 9 X 


cos 9 Z sin 9 y cos 9 X + sin 9 z sin 9 X 
sin $ z sin 9 y cos 9 X — cos 9 Z sin 9 X 
cos 9 y cos 9 x 


2 ” 

V 

Z 


frame — 1 



rame — 2 
rame — 1 


■ St ' 

Y 


( 12 ) 


( 13 ) 
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We now note three important properties of rotation matrices. 


q f rame — 2 

q/ - frame — 7\ nframe — 2 
^r/rame- 1 frame — 1 

(14) 

r nframe-lyT 
^ I'H’/ rame — i i 

r n/ram«-2iT C/ — f rame — 2 v 

= ARframe- J S ( W /rame-l) 

(15) 

SiRVrl’ZlZU) 

r> / rame — 2 ct nframe — l 

— R* frame — Nframe — 2 

(16) 


We will now derive expressions for angular velocity of the object as a function of the derivatives of the parametnza- 
tions of the orientation angles. We will now derive the denvatives of the elementary rotation matrices. 




'10 O' 

0 cos 0 1 — sin0 r 
0 sin0 x cos 0 x . 





Pi 0 0 

0 — sin0 x cos 0 x 

0 cos -sin0 T . 

'10 0 
0 cos 0* - sin 0x 

0 sin0 x cos 0x . 

= S(B s )R(9 x ) 

= S(B y )R(0y) 

= S(B t )R(B t ) 


0 0x 


0 

-0* 

0 


For the case of fixed axis rotations, we note that 




rame-2 

rame-1 


d n/rome - 2 


^ n/ram«-2 
df ^/rome-l 


d n/rome~2 

dt / ram «~ l 


— D/ rome ~^ 
^ rt /rom«-l 


S( 


-• /rom«-2 \ n/ rom«-2 
W /ram«- 1 / 

■ -/rame-2 \ 
^/rame-l ) 


5(, 


R(8,)R(e y )R{e,) 

+J R(0x)ii(tf y )[^(ex)] 

S(9.)R{9 z )R(9 y )R(9 I ) + fl(0 x )S(0 s )fl(0 y ).R(0x) 

+iJ(e,)fl(« y )S(dx)/2(9x) 

S(*,)£(0,)*(f y )*(0x) 

+R(9 z )S(9 y )R T (9 z )R(9 z )R(9 y )R(9 I ) 

+R(8 z )R(9 y )S(6 I )[R(9 z )R(9 y )\ r R(8 z )R(9 v )R(0 x ) 

[s(«x) + s(R(9 z )9 y ) + 5(i?(0z)«(« y )«»x)]^::“::! 
[S(0 X ) + S(R(9 z )9 y ) + 5(fi(0.)«((» y )<>x)]ii{"”::? 
5(9,) + S{R(9 z )9 y ) + S(rt(fi x )ii(0 y )fM 


We simplify the above expression to get, 


S(Cj 


f rom«— 2 
/ram«-l 


) 


0 — 0* o 

0, 0 0 

0 0 0 


+ 


+ 


0 0 COS0*0y 

0 0 sin0 t 0y 

— cos0,0 y — sin 0 x 0y 0 


0 

— sin 0 y 0 z 
— cos 0 y sin 0*0i 


sin 0 y 0x 

0 

cos 0 y cos 0*0 r 


cos 0 y sin0*0 x 
— cos 0 y cos 0*0x 
0 


0 

/rame-2 
x — frame — 1 
/rom«-2 
”' u ^y — /rame-1 


/ rame — 2 
z — f rame — 1 

0 

frame— 2 
^ x — J rame — 1 


frame— 2 
W y- frame — 1 
frame — 2 
~ LlJ x — frame— 1 
0 


(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 

(29) 

(30) 

(31) 

(32) 

(33) 


(34) 

(35) 

(36) 

(37) 
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frane-2 
X — frame — 1 
/rome-2 
y — frame — 1 
f rarne — 2 

z — frame — 1 _ 


cos dy cos 6 Z 9 X - sin 8 z 6 y 

= 

cos 0 y sin 9 Z B X 4 cos B z 9 y 


B z — sin B y B x 



r cos By cos B z 

— sin 8 Z 0 1 


= 

cos 8 y sin 9 Z 

cos B z 0 


— sin 8 V 


6, 


(38) 

(39) 


Note from equation (39) that determinant of the matrix relating the angular velocities of frame -2 and the derivatives 
of the parametrization cos 9 y . Therefore, the matrix is invertible for small values of the angle 8 y . 

2.4 Properties of Quaternions 

We will present here the properties of quaternion algebra that we use in this paper 1 [K.W86] [OB79], We will also 
derive the derivatives of quaternions. 

Quaternion addition: The sum of two quaternions x and y is given by 
x 4 y = 

Quaternion product: The product of two quaternions x and y is given by 


x-o + yo 
x + y 


xy = 


— X Xo 

Xol + XX X 


y = 


xoyo ~ x y 
xo y 4- yox 4 x x y 


Quaternion conjugate: The conjugate of a quaternion q is given by 

* 

q = 


9o 

~Q 


Quaternion norm: The norm of a quaternion q is defined to be 
llq|| 2 = q*q = go 2 4 g- q. 

This is analogous to the euclidean vector norm of a four dimensional vector. 
Quaternion Inverse: The inverse of a quaternion q is defined as 
1 


q" 1 - 


l|q|| 2<1 


It can be verified that this inverse has the property that 

q -1 q = qq _1 = 1 . 

Rotation operation: A rotation of a vector x by B about an axis n is given by 

qxq* 

where q is the quaternion given by 


q = 


cos(f) 

sin(|)n 


The derivatives of the Quaternion representing a rotation operation are given by 


q = q{^} 

q = q{^} + q{~^}{^} 


where w is the angular velocity and w is the angular acceleration. 

The angular velocity and angular accelerations are given in terms of the quaternions through the following relations: 
w — 2q*q 

w = 2q*q 4 2q*q 

boldface letters represent quaternions 
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Let us now calculate the derivative of the conjugate of the quaternion. As the quaternions representing a rotation 
operation are unit quaternions (unity norfi), the inverse is the conjugate. Hence 

i = qq* 

0 = qq’ + qq* 

q* = -q - 1 qq* 

= -q*qq* 

= -q*q{^}q* 

2.5 Relation between Rotation Matrices & Euler Parameters 

The operation of rotation by an angle 6 about an axis n can be represented by both a rotation matrix (R) and euler 
parameters (q) as 

R = cos(0)l + (1 - cos(0))nn T + sin(0)nx 



These two representations can be related as follows: 
R = (q 0 2 - fq)h*3 + 2 qf + 2 9o ?X 


3 . Dynamical Equations of the Rotor 

In this section we derive the equations of motion of the rigid body rotor supported by active magnetic bearings. We 
begfn thL section, with a derivation using rotation matrices, and then proceed to do the same using euler parameters 

(quaternions). 

3.1 Dynamical Equations using Rotation Matrices 

To eliminate ambiguity regarding the specification of reference frames, we will primarily work in the rotor reference 
frame, and finally transform the coordinates to the inertial reference frame. We derive the dynamic equations 
magnetic bearing in a systematic manner. For an excellent exposition on kinematics, refer to [RS94J. 

Wecompute the angular momentum of the rotor about the origin of the rotor reference frame, denoted as 8%% 


fj rotor jrotor -rotor \ / 

si rotor — * rotor rot or 

We utflize the principle of torque balance to relate the rate of change of angular momentum to the net torque. 
We note here that by net torques (T ro,or ) we refer to the summation of the applied torques (r , and the moments 
of the applied force! (f‘) about the origin of the rotor reference frame. It is to be understood of course that the 
quantities on either side of the equality will be referenced in one coordinate frame. Indeed to avoid ambiguity, we 
will henceforth refer to each of the aforementioned quantities in a single coordinate frame. Expressing all quantities 

in the inertial reference frame, we get 


rj rotor r rot or r>tnert< 

H rotor — J rotor ^rotor 


pmertial —rotor 
ttrotor ^inertial 


Step 3. 

We utilize Newton’s torque balance equations to derive 


the following. 


rjinerttat --rotor 
H-rotor ' tnertiol 



. -»rot or 
+ ^ rot or 


t=l 


X H 


rotor 

rotor 
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n 


ntnerttal — rotor 
**rotor Mnerti 


»+J2 r "° 


x R 


tnerttai 

rotor 


, ntnerttal -rotor ntnerttal nrotor Tj rotor 

“T “-rotor 'V inertial x -“-rotor ^tnerdoi^rotor 

rrotor -rotor 
* rotor ^ rotor 


, ntnertial , -rotor nrotor rrotor ntnerttal -rotor , 

' Ilr rotor l u,/ mertiad x jx xnerttai* rotor ^-rotor ^tnerttai) 
n 

ntnerttal —rotor V s C/r* \ ptncrtioi pi 

^rotor 'mertial T / ^ J 1 T rotor ) -“-rotor ‘ inertial 

tail 

4- D mert, °* Q( l ~’ i rotor \ D rotor rrotor ntnerttal -rotor 
| ■‘‘■rotor J [^inertial ) rLincrttat* rotor fl-rotor & inertial 


Trotor -rot** 
* rotor ^ rotor 


where Kotor is the point of application of the force f- nertiall and is given as 


7 X — rotor 


^x— inertial 

r 1 

' y — rotor 

_ /T-iinertial 

‘ rotor 

^y — inertial 

r z — rotor 
1 

r z — tnerttai 
1 


RV 


0 0 0 


ntnerttal 

•‘‘■rotor 


rotor — origin 
r r-»n«rtiai 
rotor— origin 
y — tnerttai 
rotor—origtn 
r z — inertioi 
1 


r i — inertial 

r y — inertial 
r » 

' z-merhal 

1 


(41) 


(42) 


... ^rotor . 

We now compute w rotor in 


the following manner utilizing (14) - (16) 


^ rotor 
^ rotor 


® r nrotor iT -rotor 
^ [■‘‘'tnertialj ^inertial 

r nrotor i T- rotor r nrotor tT a/ -rotor \-rotor 

[“•inertial J ^inertial [ “■inertial J V^tnertioi ^^incrtiai 

r nrotor -\T -rotor 
[•“•inertialj ^inertial 

nmerttai —rotor 
•‘trotor ^inertial 


Substituting 


the expression for u? rotor from equation(46) in the torque balance equation, we get 


(43) 

(44) 

(45) 

(46) 


n 

ntnerttal —rotor . V v \ nmertial ni rrotor ntnertiaf 4 rotar 

-“■rotor ‘inertial ^ J \ T rotor ) ^rotor ^inertial — ‘ rotor -t^rotor ^inertial 

isl 

i ptHcrtidi n/ -*rotor \ nrotor . rrotor ntnerttal —rotor 
i •‘‘•rotor *- 3 \ t4, 'inertiat /l-“ , mertial-‘rotor •‘‘■rotor ^inertial 

We now recast the above equation in the following form 


n 


^rotor 

^inertial 

rrrotor n»n«rtio(i — l r ninertiol -rotor i \ ^ Of 2 4 \ ntnerttal i 

[ J rotor -“-rotor J i^rotor 'tnertiaJ ' / t -- 1 1 ‘rotor J^rotor -‘mertiaij 

(47) 


\ jrot or ntnerttal} — \r ntnerttal q/ -rotor \ nrotor rrotor ntnerttal -rotor i 

L- 1 rotor ■‘‘■rotor J [•‘‘■rotor v* inertial ) **’\nert\al* rotor -^rotor ^inerttai J 

(48) 

Rrotor 

^inertial 

__ nrotor r rrotori-l r ntnerttal -rotor , \ ' c / \ ninertiai ni •» 

— “-inertial L-* rotor J [‘Votor “inertial T / ^ [‘rotor jH-rator -‘inertialj 

(49) 


nrotor rrrotori-1 nmertial c/.-l^or ^nrotor rrotor ‘ntnerttal -rotor 

■“■inertial [-* rotor J rotor J \ Uf inerttal ) t*-tnertial*rotor -K-rotor ^ inertial 

(50) 

Step 4. 




We derive the force balance equations by first calculating the expression for the linear momentum of the rotor 
in the following manner. 


Z rotor rotor —rotor 

rotor — 771 v rot or 

rotor ntnerttal —rotor 

TO IX rotor ^»n«rtiol 

Step 5 

Newton's law asserts that the rate of change of linear momentum equals the applied force. That is, 


n 

(inertial \ ' pi 

■rotor / j -^inertiai 

i=l 

^ rfrotor-i 

— ^L^rotorj 

(51) 


rotor - rotor - rotor Trotor 

TO v rotor + ^ rotor x — rotor 

(52) 


rotor - TOtOT n/ ntnerttal —rotor \ rotor ntnerttal -rotor 

— TO Nrotor i ■‘trotor ''•‘-'inertial **■ rotor ^inertial 

(53) 
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We now compute vlTtZ in (53) in the following manner utilizing (14) - (16). 


jrotor 
^rotor = 

d e nrotor tT -rotor 

TTl-tt, nertial] u mertial 

at 

(54) 


rnrotor , T - rotor rnrolor 1 T C(^ rotor . 

[/t mer tialj ^inertial nertial] ^ t ^inertial / v mertial 

(55) 

= 

ninertiai - rotor ninertial q{ -rotor x -rotor 

“rotor ^inertial ^rotor * “^inertia! / '•'inertial 

(56) 

Substituting the expression for v7JtZ from equation (56) in the force balance equation (53), 

we recast the force 

balance equation 

as, 


n 

. ; V — > -t. rotor ninertial - rotor rotor ninertial r v -rotor \ -rotor 

R'7oYJr ^ ^inertial = m ^rotor inertial “ m ^rotor *S(w inertia , )*, nert.ai 

(57) 


rotor c / ninertial -.rotor x ni nertial -rotor 
-f-m O ^ t*-T 0 tor ^inertial f^rotor ^inertial 


(58) 


We now rewrite the force equation (58) in the following form 


•jrotor 

finertiai 


n 

1 V - '' ft* r./ -.rotor -rotor 

ro "— > rin«rtial ' •3\ w mert*ai/ w mertial 


nrotor c / ninertial -rotor \ ninertial -rotor 
— ^tnertiat^l^rotor ^tnertiaU^^otor ^mertiai 


Note that the last term in equation (61) may be simplified in the following manner, 


nrotor 

^mertioi 


C/ ninertial -rotor \ ninertial —rotor 
•^V-ttrotor ^ inertial ) ^rotor u inertial 




rotor sr nrotor ]T 

rotor /l^mertioi j u mertia< 
— rotor \ —rotor 


c/ nrotor -rotor \ -rotor 
O \ itmertiol^rotor j u inertia< 
Qf —rotor \ — rotor 

*■> (^inertial /^inertial 


Substituting the expression in equation (64) in equation (61), we arrive at the force balance equations, 


nrotor 

u m«rt«ot 


E f<» , n/ -rotor \ -rotor c/ -rotor \ -rotor 

F t nertial + *5 (^merttal Mnertiai “ l w mertiaf /^inertial 


7 7f57 £ ^' ntr " ai 


(60) 

(61) 


(62) 

(63) 

(64) 


(65) 

( 66 ) 


Collecting 

the force and torque balance equations, we write the dynamic equations of the magnetic bearing as 

follows. 



jrotor 

^inertial 

n 

nrotor r jrotor'i — X r ninertial ^rotor , \ O/J 4 ^ ninertial pi 1 

= RinertialUrotor] [Rrotor r mertial + / , J ( r rotor J ^rotor ^mertialj 

(67) 


nrotor rrrotori -1 ntnertial c/ -rotor x nrotor rrotor ninertial -rotor 

'*^inertial[frotorj ^rotor •jfWmertial )dv t ner*ial4 rotor ^rotor inertial 

(68) 

■jrotor 

^inertial 

- m ro<or £ ^inertial 

1=1 

(69) 


3.2 Dynamical Equations using Quaternions 

As we saw in the previous section in equation 69 the force balance equations essentially are a restatement of F = ma 
in the inertial coordinates. So we will only consider the angular momentum equations. 


rprotor 
* rotor 


n 


-rotor 
T rot or 


+ 


£ 


r rotor 


i=l 


X P* 

A 1 rotor 


d r rr rOtOr 


} 


T rotor -rotor -rotor 
* rotor ^ rotor ■ ^ro tor 


x H\ 


rotor 

rotor 


Let us look at the first term on the right hand side of 71 initially. 


jrotor nrotor 
*rotor ^ rotor 


.-rotor d . inertial -rotor inertial * 
i rotor TT (Qrotor w tnerttaiMrotor 

at 


rrotor / « inertial —rotor inertial* 
i rotor (Qrotor W, nertia <q r otor 


inertial ~ rotor inertial * 
T Qrotor ^inertialMrotor 


( 70 ) 


344 


i inertial -rotor -inertia/* 
"T* Hrotor “■'inertia/Hrotor 


jTotor / inert! at t 1 -rotor | -.rotor inertia/* 
1 rot or t Hrotor \ ^ ^inertial J ^inertia/ Hrotor 


. inertial p rotor inertial* 
* Hrotor ^mertta/Qrotor 


inertial -»rotor inertial* inertial t 1 -*rotor \ mertia/*\ 

Hrotor ^mertia/Hrotor Hrotor 1 ^ ^tnert la/j Hr otor ) 

_ rrotor/ inertial rrotor -i -# rotor inertial* , inertial ^ rotor inertial* 
J rotor\Hrotor \ ^ w inertial J^inertial Hrotor ' Hr otor ^inertia/Hrotor 


inertial -rrotor r^-«rotor i inertia/* 

Hrotor ^inertia/ l ^ ^inertial /Hrotor 

rrotor / inertial r*rotor inertial* \ 

1 rotor vHr otor ^mertia/Hrotor ) 


Substituting this back into 71 we get 


n 

-•rotor . \ ^ -4' 

'rotor 1 y ^rotor 

i=l 


Pro 


rrotor / inertial rrotor inertia/ 
J rotor \Hrotor ^inerti a/Hrotor 


) -f W 


rotor 

rotor 


X 


a 


rotor 

rotor 


(71) 


inertial -trotor inertia/* , \ t 

Hrotor ^inertia/ Hrotor T ) ^ (Hr 


rotor ^inertia/ Hrotor 


inertia/ pi 
A Hrotor 1 u 


inert«aZ*\ 
nertia/Hrotor ) 


rrotor / inertial m ? r °t° r inertial* 

* rotor \Hrotor ^ inertia/ Hrotor 


, inertia/ -jrotor inertia/* inertial Arotor inertia/* 

■ Hrotor ^mertia/Hrotor Hrotor inertia/Hrotor 


(72) 


7*rotor 

^inertia/ 


inertia/ , *r rrotori — 1 r inertial -jrotor V^Y_ inertial inertial* _ inertia/ £S» 

| TOIOT [ i rotorJ lHrotor r inerttaf ~T / ^ Hrotor r t nertia/Hrotor * Hrotor ^inertia 


_ inertial -»rotor inertial* inertial itrotor -i 

Hrotor ^mertia/Hrotor ^ Hrotor ^tnertta/J 


(73) 


Comparing equations 68 and 73 we see the equivalence of both these derivations. 

We have briefly derived the dynamics equations of the magnetic bearing system using euler parameters, as they 
have many advantages. Euler Parameters are defined everywhere and they have a nonsingular mapping with the 
rotational velocity w . Using Quaternion algebra the above expressions can be further simplified. Simple expressions 
for all composite rotations and rotating reference frames can be developed [K.W86]. Euler parameters have also been 
shown to be as efficient computationally as rotation matrices and more compact in storage [JR90]. 


4 Small Angle Assumption 

We have derived a detailed nonlinear model of the rotor supported by active magnetic bearings. We will now present 
the standard assumptions made in deriving the dynamical equations of the magnetic bearing and the simplification 
achieved on the nonlinear model [FK90]. 

In the magnetic bearing system, let the spin axis be x and the pitch and yaw axes be y and z axes. Let the 
spin angle, pitch and yaw angles be 9 x ,9 yy 0 x . Usually we assume that the angles 9 y ,9 z are very small so that 
cos(^),cos(^ z ) % 1, sin(#y), sin(0 z ) 0. Also we can reasonably assume that the product of velocities and angular 
velocities are small and can be ignored. The external forces acting on the system are the forces at the two radial 
bearing systems F y , F z and F yi F z \ the force at the axial bearing F z and the external torque r mo /or applied along 
the spin axis. With these assumptions the equations of motion of a rotor supported by magnetic bearings reduce to 


X 


X 


y 


y 


z 


z 


9x 


9x 


F x 

9y 


0y 


C 1 

4 y 

9z 

X 

= F 

9z 

X 

+ G 

Fl 

y 


y 


Ft 

z 


z 


Tmotor 

9x 


9x 


Py 


9y 


. 9z _ 


9z 
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where F £ 3ft 12x12 and G £ 3£ 12 * 6 . 

06x6 ^6x6 

0 0 0 0 0 0 

0 0 0 0 0 0 

F = 00000 0 

° 6 * 6 0 0 0 0 0 0 

0 0 0 0 0 -a>a 

0 0 0 0 ua 0 

where a 6 is a constant dependent on the inertias J*, J y 


Feedback control of the rotor motion 


The rotor system motion is decoupled between the spin axes and its pitch and yaw axes.Hence for the design of a 
linear state feedback controller we shall consider the dynamical equations of only the pitch and yaw axes motions of 
the rotor system, given by 


04x4 ^4x4 

04x4 -A[u/] 


where x 


[x\x\x\x\x\x\x%x\) T G » 8 , « = [ttVuVf G 3? 4 , A : 


4 , B G » 4 


= 


0 0 0 0 

0 0 0 0 

000 -wfl 

0 0 wa 0 


■h -h 0 0 

0 0 MM 

0 0 —j j 

Jy Jy 

1l -±l o 0 


Note that the x\,x\ subsystem and x\,x\ subsystem form two decoupled 2x2 systems, while x 3 u x \,x\,x* 2 form 
a coupled 4x4 system. Let us choose the control 
f -k\x\-k\x\ ' 

i —kjif — k\x\ (76) 


u = B 


— kl Xi — klx 2 

*S*a 

-k\x\-kUi 


such that k 1 > 0 fori — 1,2 : j = 1 • • *4. We can now show that this control stabilizes the x\ , x-t ,x\,x\ system at all 
speeds u>. 

Theorem 5.1 

Given (Gl) The system given by equation 1 4 
(G2) Feedbacks u given by equation 76 
Then (Tl) u\u 2 stabilize the system at all speeds u>. 

Proof: 4M> Let us choose a Lyapunov function candidate V as follows: 

Ed 

V 2 

Taking the derivative of the Lyapunov function we get 


V = F, k l X l X 2 ~ * 1 * 1*2 - kjx 3 2 


3 4 3 4 

— (jjax 2 X2 + u)ax 2 x<2 


= ^-kfri 


< 0 ifx 3 2 fcOforj = 1 •••4 
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Figure 1: y vs time : w = 0 



Figure 2; O y vs time : w = 0 


Hence by Lyapunov theorem x> 2 converge to 0 for j = 1 • ■ - 4. As the maximum invariant set containing the set x{ = 0 
is x 1 - 0, by LaSalle theorem [M.V78] x] also converges to 0. Hence the system is stable at all u. 0 4 

5.1 Simulation Results 

We simulate this system using the control system simulation package SIM NON. We present the simulation results of 
applying the state feedback control given by equation 76 to the magnetic bearing system. We simulate the system 
response y and 0 y under the feedback when there is no spin (Figures 1,2) and with a spin of lOOrad/sec (Figures 3, 



Figure 3: y vs time : w = 100 
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5.2 Magnetic Bearing Model 

In this section we present the basic equations for a single axis magnetic bearing, and two associated pole excitation 
schemes. The dynamical equations of the magnetic bearing may be written as follows. 

(77) 

4>o = « 0 (7g) 


<Pn 

= XL 

tr 

4>0 


flux at pole location 0 

4>n 

G » 

flux at pole location n 

UQ 

€ » 

Control at pole location 0 

U n 

G R 

Control at pole location 7r 


Let the control force F generated by the magnetic bearing be the net force produced by the bearing elements at 
the angles 0 and it (the positive and negative poles in a pair). Indeed, 

(83) 

F = Fo - F„ 

We shall design the feedback control of the rotor using the net force as the control actuation. Treating th,s 
requisite force as the commanded output of the magnetic bearing subsystem described by equations 77, 78, we design 
the flux feedback as a deadbeat controller. Inherent in this approach is the assumption that the flux feedback loop 
would be run at a much faster rate than the bandwidth of the force feedback system. 

Discretizing the flux equations in the following manner 

*,(* + !) 1 _ M*) 1 . [ Tuo(k) 1 (84) 

<Mfc + !) J L <M*) J l Tu *( k ) . 

where 

(85) 

un(k) = is the net control voltage at pole 0 

1 ^ , , ( 86 ) 
u n (k) = is the net control voltage at pole tt 

We now note the relation between force and flux is given the following form 

T(k+ 1) = A'/orce-/lui[0o(* + 1) “ <t>l( k + J )] ^ 

where the magnetic constant « 6 » elating the forces produced due to fluxes applied is assumed to be 

known. Choose the control inputs in equation (84) to be of the following form. 

—<t>o(k) + vo (k) 

tio(k) - j. 

- 0 r (k) + v w {k) 
u„(k) = j 
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( 88 ) 

(89) 



where 


vo(k) € K is an exogenous control input, specified later 

r M*) € ft is an exogenous control input, specified later 

Substituting equations (88) and (89) in (84), we get, 


<z>o { k + 1 ) 


■ vo(k) ' 

Ow( k -f* 1 ) 


[ ».(*) . 


Substituting equation (92) in (87), we get, 

F(k + 1) = K force- flux[vl(k) - ?4(*)] 

We consider the following choices for choosing the control inputs v 0 (k) and v n (k). 


(90) 

(91) 


(92) 


(93) 


5.2.1 Mutually Exclusive Scheme 
Choose the control to be 


w(k) = 


~F(k + 1) 

A f or ce- flux 


Now choose the controls v 0 (*) and v w (k) in the following manner. 
w(k)> 0 - / Mk) = \F4k) 

’ { M*) = o 


w(k) < 0 


wo(t) = 0 
v„(k) = ,/|u>(*)| 


(94) 


(95) 

(96) 


5.2.2 Biasing Scheme 

Choose the following structure for the controls v 0 (k) and v*(k). 

Uo(l“) — H" Vvana6fe(£) 

V*{k) = Vvano6le(l) 

where 


(97) 

(98) 


€ 3? is a constant biasing input 
Vvanabie € 5? is a varying control input 
Note that such a structure for vo(k) and Vrr(k) implies that 

A/orce-/lu*[vo(*) “ *£(*)] = K f orce- f tux[{v bia . + V vartabte ( *) ) 2 - [(vfei a , - V W aria«e( *))*] 

^ / orce — f lux[^VbiasVvartablek] 

We now choose the control v va r,abu(k) to be 


(99) 

( 100 ) 


( 101 ) 


Vt»orio6i«(^) — 


- gCj+j) 


K 


force — flux Vfci at 


( 102 ) 


where the control F is chosen to stabilize the rotor motion. 

Both these excitation schemes have their advantages and disadvantages. In the constant biasing scheme, we note 
that the force to flux equations become linear. Also by choosing as the control is scaled by v blaa , change in force (or 
equivalently currents) required for a certain net force can be reduced. But maintaining a constant biasing voltage may 
increase the losses. An alternative might be to use permanent magnets to provide the bias voltage. In the mutually 
exclusive scheme we provide a force (or current) in only one pole, from a pair, at any time. On the other hand, the 
force to flux relations are nonlinear. 


6 Multiple Sensors &: Redundant Actuators 

In many situations, we measure the same output with multiple sensors and the measurements have to be averaged in 
some manner. Similarly, in the case when we have redundant actuators (more than the necessary three orthogonal 
pairs), we need to apportion the actuation forces in an optimal sense, between all the actuators. Linear least squares 
theory provides us with a method for doing these [RH88] [Aub79] [J.L55], In this section we will look at using the 
least squares estimation schemes for averaging measurements from multiple sensors and splitting the forces among 
redundant actuators. 
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6.1 Linear Least Squares 

Definition 6.1 .4 complete inner product space X is called a Hilbert space. 

Definition 6.2 Given a Hilbert space X, and a subset U € X, the orthogonal complement of U, denoted by U L is 
defined aj follows . 

{/■*■ = {x € X .< x.u >= OV u € U) ,103) 

That is. the orthogonal complement of a set U € -Y is the set of all vectors in -Y that are orthogonal to every vector 
in U. 

Theorem 6.1 Projection Theorem 
Given (Gl) .4 Hilbert space X . 

If (It) J £ X is a closed subspace of X . 

Then (Tl ) The Hilbert space X can be decomposed into the direct sum, 

X = U®U L (104) 

Definition 6.3 Let U 6 X be a closed subspace of a Hilbert space X. Decompose a vector x € X into the direct sum 
x = fo + fi where x 0 € U and x x € . Then x Q is called the orthogonal projection of the x $X onto the subspace 

U 6 X . 

Theorem 6.2 Projection Theorem 
Given (Gl) A Hilbert space X. 

G(2) A direct sum decomposition of X — U © U x . 

G(3) .4 vector x € X . 

If (II ) xq is the orthogonal projection of x onto the closed subspace U € X. 

Then (Tl ) x — xo «J the orthogonal projection of x onto the close subspace U 

Theorem 6.3 Minimum Norm 
Given (Gl) .4 Hilbert space X, and a vector x € X. 

(G2) A closed subspace U 6 X. 

If (II) fo is the orthogonal projection of x onto the subspace U . 

Then (Tl ) For each u 6 U , 

||x - roll < III - “II ^ 105 * 

Given two Hilbert spaces X, Y, let the operator A be such that A : X — • Y. We now make the following 
definitions. 

Definition 6.4 The range of At X — Y denoted as U(A) = {/l[r] € Y VrfX). 

Note that the range of A is the set of all vectors in Y that are obtamed by the action of the operator A on every 
element in X , That is t TZ-(A) C Y . 

Definition 6.5 The null space of A : X — Y denoted jV(A) — {x g X : ^l[x] = Oy } 

Note that the null space of A is the set of all vectors in X that are mapped by A into the zero element of Y. It is 
clear that A/*(x4) C X. 

Definition 6.6 The adjoint of a linear operator A : X Y , denoted as A*, is defined as follows. 

• A* : Y X 

• < A[x], y >y = < x,j4*[y] >x V x € X, y € Y. 

where < >x *s the inner product defined in space X, and < >v it the inner product defined in space Y. 

The usefulness of the adjoint operator will become evident in the solution of Unear equations. The foUowing properties 
of the adjoint operator are vital to its use. 

. Given an operator on a Hilbert space A : X - Y, and its adjoint A* : Y - X, it can be shown that 

1. Af(A) = M(A*A) 

2. H(A) = Tl(AA*) 

• Given an operator on a Hilbert space A: X — Y, and its adjoint A* :Y — X, it can be shown that there exist 
orthogonal direct sum decompositions of Hilbert spaces X and Y of the foUowing form. 

1. X = K(A*)®Af(A) 

2. Y = K{A)®Af(A*) 
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6.2 Least Squares Solution of y — A[x] 

Given a linear operator on the Hilbert space A : X — ► Y, and a specific y\ £ Y, we define the solution of the linear 
equation y\ — A[x] as {x £ X : y\ = A[x]}. 

There are three cases that merit consideration. 

• If the operator A : X — ► Y is such that the 71(A) = T and the Af(A) = {Ox}, then the solution of yi = A[x] 
exists and is unique. The solution is given as x = A~ l [y\\. Note that the inverse A~ 1 : Y — ► X exists. Such a 
solution corresponds to a system of linear equations with as many equations as there are unknowns. 

• If the operator A : X — ► Y is such that the 71(A) C Y and the A f(A) = {Ox}, then we note the following, 

yi = A[x] 

A*[y i] = A*[A[x]] 

A*[ yi] = A*A[x] 

where the operator A* A : X — ► X. Note that Af(A*A) — Af(A) = {Ox}. This implies that the inverse 
(A*A)“ 1 : X — ► X exists. The solution therefore can be written as 

x =r (A'Ar'A'yt (106) 


There is a simple geometric interpretation of the above result. Given y\ 6 Y } there is a unique direct sum 
decomposition of y\ as, y\ = (yn 6 71 (A)) 0 (y ^ 2 € Af(A*)). That is, the vector in 71 (A) closest to y is the 
vector y — 3/12. Indeed, the best one could do is to find a solution x € X such that A[x] = y — y 12 = y 11. So 
we attempt the following solution, 


y\ - Vi2 = A[x] 

A*[yi - yi 2] = A* A[x] 

- A*[y\ 2 ) — A*A[x] 

,4*[t/i] — 0y = A*A[x] 

A*[yi] = A*A[x] 

x = (A*A)- l A*{ yi ] 

The solution (106) is called the least-squares solution of the linear equation y\ 
corresponds to an overdetermined set of linear equations. 

• Given a linear operator A : X — * Y is such that the 71(A) = Y and the {Ox} C Af(A), we follow the geometric 
intuition as follows. 


A[4 Such a solution 


- Solutions exist as 71 (A) = Y . 

— Consider any solution Xi € X : y\ = v4[z;]. This solution has a unique direct sum decomposition of the 
form x t = (x 1 1 6 7 l(A*)) 0 (x,2 € A/^(A)). Indeed, there is no contribution of x,2 € A 7 (A) to the solution 
of yi = A[x). Furthermore, as x,i € 7t(A k ) i it is true that there exists w £ Y such that x,i = 

Note that 


y 1 = A(x x ] (107) 

= A[xn + x t2 ] ( los ) 

= A[xn] (109) 

= A[A*[w)] (110) 

= AA*[w] (111) 


Now note that AA * : Y — ► Y. Also 7 l(AA*) = 71 (A) = Y. This implies that Af(AA*) = 0y. This 
guarantees that (AA*) -1 exists. We therefore solve for w in equation (111) as 


w = (A A*) 


( 112 ) 


Note that the minimum norm solution is certainly one that does not include elements from .V*(A). There- 
fore, the minimum norm solution of yi = A[x] is x,i = A*[w] = A*(AA*) -l j/i . This solution corresponds 
to an underdetermined set of linear equations. 


6.3 Least squares solution to multiple sensors and redundant actuators 

Let us consider the case when there exists a multiplicity of sensors for the same measurement. Let the actual 
measurement we are looking for be x and the multiple sensor measurements be y = Ax. Then, to get a mean 
measurement, with minimum error to the actual measurement, corresponds to exactly the overdetermined case in the 
least squares estimation. The measurement is then given by 
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This (A*A)^ 1 A* is indeed the pseudoinverse of the A matrix. 

Now consider the case when we have redundant sensors and we are looking for a force split that minimizes the 
norm of the total force applied. Given forces x produced by the redundant actuators, the net force applied is given 
by y = Bx. Now, given a force to be applied y , splitting it among the redundant actuators with minimum norm, is 
exactly the underdetermined case derived in the least squares estimation. The solution is given by 

* = B*[w] = B*(BB*)~'y 

Note that is the pseudoinverse of B . 

7 Summary 

In this paper we have developed the detailed dynamical equations of a rigid body rotor supported by actively 
controlled magnetic bearings. We have done this using both Rotation Matrices and Quaternions to see the equivalence. 
Quaternions are more convenient to use, as they provide a nonsingular (invertible) transformation to the angular 
velocity w. Also euler parameters are computationally as efficient and more compact in storage than rotation matrices. 
In addition, in developing the model of the magnetic bearing system, we have oonsidered two schemes for pole 
excitation. 

We notice that the model of the bearing system depends on the angular velocity in the spin direction. We have 
developed a state feedback controller that stabilizes the system for ail speeds of rotation. We also note that this 
controller essentially decouples the system into 2x2 subsystems. We have presented simulation results showing the 
performance of the controller. 

Finally we also present a least squares scheme for minimizing the residual in measurements of output with multiple 
sensors, and for minimizing the norm of the actuation forces when there are redundant actuators. 
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